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Abstract 

The Calogero model bears, in the continuum limit, collective excitations in the 
form of density waves and solitary modulations of the density of particles. This sector of 
the spectrum of the model was investigated, mostly within the framework of collective 
field theory, by several authors, over the past fifteen years or so. In this work we shall 
concentrate on periodic solutions of the collective BPS-equation (also known as "finite 
amplitude density waves"), as well as on periodic solutions of the full static variational 
equations which vanish periodically (also known as "large amplitude density waves"). 
While these solutions are not new, we feel that our analysis and presentation add to the 
existing literature, as we explain in the text. In addition, we show that these solutions 
also occur in a certain two-family generalization of the Calogero model, at special points 
in parameter space. A compendium of useful identities associated with Hilbert trans- 
forms, including our own proofs of these identities, appears in Appendix A. In Appendix 
B we also elucidate in the present paper some fine points having to do with manipulating 
Hilbert transforms, which appear ubiquitously in the collective field formalism. Finally, 
in order to make this paper self-contained, we briefly summarize in Appendix C basic 
facts about the collective field formulation of the Calogero model. 
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1 Introduction 



The Calogero Model (CM) [1] - [3] is a well-known exactly solvable many-body 
system, both at the classical and quantum levels. It describes particles (consid- 
ered indistinguishable at the quantum level) on the line, which interact through an 
inverse-square two-body interaction. Its quantum Hamiltonian is 

H_ If I A(A-l) f 1 

where m is the particles' mass, and the dimensionless coupling constant A parametrizes 
the inverse-square interaction between pairs of particles. ^ 

The CM and its various descendants continue to draw considerable interest due 
to their many diverse physical applications. A partial list of these applications can 
be found, for example, in the introductory section of [4]. For recent reviews on the 
Calogero- and related models see, e.g., [5, 6]. In addition, for a recent review on 
the collective-field and other continuum approaches to the spin-Calogero-Sutherland 
model, see [7]. 

In the present paper we concentrate on the thermodynamic limit of the CM. In 
this limit the system is amenable to large- TV collective- field formulation [8, 9, 10]. 
As is well-known, the collective theory offers a continuum field-theoretic framework 
for studying interesting aspects of many-particle systems. Clearly, a description of 
the particle systems in terms of continuous fields becomes an effectively good one in 

the high density limit. In this limit the mean interparticle distance is much smaller 

^Note that we did not include in (1) a confining potential. This is not really a problem, as we 
can always add a very shallow confining potential to regulate the problem (in the case of piirely 
repulsive interactions), or else, consider the particles confined to a very large circle (i.e., consider 
(1) as the large radius limit of the Calogero-Sutherland model [2]). We shall henceforth tacitly 
assume that the system is thus properly regularized at large distances. 
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than any relevant physical length-scale, and the 5-function spikes in the density field 
(4) below can be smoothed-out into a well-behaved countinuum field. All this is in 
direct analogy to the hydrodynamical effective description of fluids, which replaces 
the microscopic atomistic formulation. Of course, the large density limit means that 
we have taken the large- N limit, as was mentioned above. 

The collective-field Hamiltonian for the CM (1) is given by [11] 

Hcoii = ^ / ^x7r(x) p{x) ^^^(^) + ^ / P(^) {~2~~p~ ^i~x^^y) "^^^'"^ ' 

(2) 

where Hsing denotes a singular contribution [12] 



HM„n = jdxp{x)d, ^ 



2m J X — y 



and P is the principal part symbol. 
Here, 



4m 

y=x 



dxd^ 5{x-y)\ , (3) 



N 



p{x) = Y.5{x-Xi) (4) 



i=l 



is the collective - or density - field, and 

6p{x) 

is its canonically conjugate momentum. It follows trivially from (4) that the collec- 
tive field is a positive operator 

p(x)>0, (6) 
and that it obeys the normalization condition 

oo 

j dxp{x) ^N. (7) 
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The latter constraint is implemented by adding to (2) a term /jL yj dx p{x) — Nj , 
where is a Lagrange multiplier (the chemical potential). 

The first term in (3) is proportional to p{x). Therefore, its singular coefficient 
— ■^dx amounts to a shift of the chemical potential a by an infinite constant. 

The last term in (3) is, of course, a field independent constant - an infinite shift of 
energy. In order for this paper to be self-contained, we have briefly summarized the 
derivation of the collective-field Hamiltonian (2) in Appendix C. 

It is worth mentioning at this point that the Calogero model enjoys a strong- 
weak-coupling duality symmetry [13, 14]. At the level of the collective Hamiltonian 
(2), these duality transformations read 

A ' ^ ~X ' " ~ A ' ^^^^ " , and 7r{x) = — , (8) 

and it is straightforward to see that these transformations leave (2) (including the 
chemical potential term) invariant. The minus signs which occur in (8) are all 
important: We interpret all negative values of the parameters and densities as those 
pertaining to holes, or antiparticles. Thus, the duality transformations (8) exchange 
particles and antiparticles. (For more details see e.g. Section 3 of [4], and references 
therein.) 

It is well-known [9] that to leading order in the expansion, collective dynamics 
of our system is determined by the classical equations of motion resulting from (2). 

The simplest solution of these equations is the constant condensate p{x) = pQ (and 
7r{x) = 0) corresponding to the ground state. 

More interesting solutions of these equations include various types of periodic 
density waves and soliton configurations [15, 16, 17]. As we explain in Section 2 



below, these periodic density waves can be thought of as a crystal made of the 
localized soliton solution. 

Recently, density wave configurations of this type were studied, among other 
things, in [18], where a certain regulator, first introduced in [16], was used to tame 
the effective collective potential. The BPS-equations associated with the regulated 
potential were then converted into a Riccati equation, which was solved explicitly. 

Such static periodic density waves are the focus of the present paper as well. 
As in [18], we convert the BPS-equations associated with the equations of motion 
of (2) into an explicitly solvable Riccati equation. However, unlike [18], we avoid 
introducing any unconventional regulators in (2). In addition, we also construct non- 
BPS solutions of the equations of motion, which are simple shifts of BPS-solutions 
by a constant, and compute their energy densities. That constant is fixed by the 
equation of motion and turns out to be either the maximum or the minimum value 
of the corresponding BPS-solution. Thus, these non-BPS density profiles vanish 
periodically and coincide with the large-amplitude waves reported in [17], albeit 
without too many details of their construction. We beheve that the constructive way 
in which we derive our static periodic BPS and non-BPS configurations complements 
the discussion in [15, 17, 18]. Since these non-BPS solutions vanish periodically, we 
can also refer to them as vortex crystals, as they constitute a periodic generalization 
of the vortex solution of [16]. 

In the present paper we also show how these known solitary and periodic wave- 
solutions appear in the collective field theory of the two-family generalization of 
the CM, under very special conditions on the coupling constants. The two-family 
Calogero model is a generalization of (1) into two species of identical particles. The 
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Hamiltonian of this model reads [19] 



1 ^ A2(A2 - 1) ^ 



2^2 ^1 ax„2 2m2 ^ (xc, - a;^) 

1/1 1 \ iVl iV2 -1 

+ 2tH--)MA.-l)gi:(j-^. (9) 

Here, the first family contains A^^i particles of mass mi at positions Xi, i ~ 
l,2,...,Ni, and the second one contains A^2 particles of mass m2 at positions 
Xa, OL = 1,2,...,A^2- All particles interact via two-body inverse-square potentials. 
The interaction strengths within each family are parametrized by the coupling con- 
stants Ai and A2, respectively. The interaction strength between particles of the 
first and the second family is parametrized by A12. 

In (9) we imposed the restriction that there be no three-body interactions, which 

requires [19]- [24] 

Ai A2 A12 ^^Q^ 



mi^ m-^ m\mi 
It follows from (10) that 

A?2 = A1A2 . (11) 
We assume that (10) and (11) hold throughout this paper wherever we discuss the 
two-family CM. The Hamiltonian (9) describes the simplest multi-species Calogero 
model for particles on the line, interacting only with two-body potentials. 

In [4] we studied the collective field theory of the two-family CM. The corre- 
sponding collective Hamiltonian is 

Hcoll = f dx d^7^i{x) pi{x) d^7ii{x) 



+ 



2mi J \2 pi J X — y J x — y ^ 

— / (ia;5a;7r2(x) p2(a:;) 9j;7r2(x) 

1712 J 

2m2 J p2 J X — y J X — y 



+ Hsing , (12) 

which is a straightforward generahzation of (2). Here Pa{x) are the collective density 
fields of the ath family (a = 1,2), and 7ra{x) are their conjugate momenta. As in 
(2), the term Hsing denotes a singular contribution which is a straightforward 
generalization of the one- family expression (3). Given that there are Na particles in 
the ath family, the densities Pa{x) must be normalized according to 

oo oo 

j dx pi{x) = Ni , j dxp2{x) = N2 ■ (13) 
—00 —00 

As in the one-family case, these normalization conditions are implemented by adding 

/ 00 \ 
to (12) the chemical-potential terms Z)a=i,2/^a I / dx pa{x) — Na \ ■ 

As was discussed in [4], the collective Hamiltonian (12) is invariant under an 
Abelian group of strong-weak-coupling dualities, which is a generalization of the 
single- family case (8). A remarkable consequence of these duality symmetries (see 
Section 3.1 of [4] for more details) is that when one sets 

A1A2 = 1 , A12 = -1 (14) 

in (10), the two-family CM (12) becomes similar, in some sense, at the level of col- 
lective field theory, to the original single family CM, with a collective Hamiltonian 



effectively given by (2), for a single effective density. More precisely, this similarity 
manifests itself in the fact that at the special point (14), the original Hamilto- 
nian (12) can be mapped by these duality symmetries onto a two-family collective 
Hamiltonian in which the two families are still distinct, but have common mass and 
two-body interaction couplings, and therefore, the two densities can be combined 
into a certain effective one-family density Pe//- In fact, at these special points, the 
classical densities (i.e., the static solutions pi{x) and p2{x) of the equations of mo- 
tion associated with (12)) turn out to be proportional to each other, and of opposite 
signs. Thus, for example, for 7712 = < 0, the common parameters mentioned 
above are A = Ai and m = mi, leading to an effective one-family density 

Peff = Pi - , (15) 

which satisfies the static equation of motion of the single-family model (2) with these 
common parameters, whereas for mi = — ^ < 0, one obtains similar relations, but 
with the two families interchanged. (Negative masses and densities in these formulas 
are interpreted as quantities corresponding to holes rather than particles, as was 
mentioned above.) 

In conclusion of this introduction it is proper to mention that the Heisenberg 
equations of motion of the collective field p{x) and its conjugate momentum 7r{x) 
may be interpreted as the isentropic hydrodynamic flow equations of an Eulerian 
fluid [25] (see also [10] ) and the latter may be associated with the completely 
integrable and soliton-bearing Benjamin-Ono equation, both at the classical level 
[26, 27] and the quantum level [28] . 

This paper is organized as follows: In Section 2 we solve the static BPS equation 
associated with the one-family collective Hamiltonian (2) by converting it into a 
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Riccati equation which wc then solve exphcitly. The solution is a static periodic 
density wave - the finite amplitude wave solution of [15]. Conversion of the BPS 
equation into a Riccati equation is achieved by considering the complex valued 
resolvent ^{z) associated with the positions of the N particles on the line (see 
Eq.(20)), whose boundary value, as the complex variable z approaches the real axis, 
is a linear combination of the density field p{x) and its Hilbert-transform p^{x) (see 
Eq. (22)) [18, 28]. That the latter combination satisfies the Riccati equation then 
follows from the BPS-equation and its Hilbert-transform. We then study various 
limits of the periodic solution. We conclude Section 2 by showing that the coupled 
BPS-equations, associated with (12) at the special point (14) in parameter space do 
indeed collapse into a single-family BPS equation. 

In Section 3 we consider the static limit of the equation of motion associated with 
(2) - namely, the full variational equation. Every solution of the BPS-equation is, of 
course, a solution of the full variational equation. It is more challenging to find non- 
BPS solutions of the latter. We seek such solutions in the form of BPS configurations 
shifted by a constant, as was mentioned above. For each of the cases A > 1 and 
< A < 1 we find two types of solutions, namely, a positive periodic density wave 
(a vortex crystal) and a negative one (an anti- vortex crystal). We discuss how these 
solutions map onto each other by the duality transformations (8). Then, we discuss 
the energy density of these non-BPS solutions, averaged over a period. We end 
Section 3 by showing that the coupled variational equations, associated with (12) 
at the special point (14) in parameter space collapse into a single- family variational 
equation. 

For the sake of completeness, and also for future use, we provide and prove in 
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Appendix A a compendium of useful identities involving Hilbert-transforms. 

In Appendix B we note and also resolve a mathematical paradox associated with 
the variational equation. It has to do with the trilocal term in the density fields ob- 
tained by expanding the square in (2). In many papers on the collective approach to 
the Calogero model, that trilocal term is converted into a local p^{x) term by employ- 
ing a certain identity among distributions, Eq.(A.20). However, strictly speaking, 
that identity is valid only for distributions acting on test functions which are inte- 
grable along the whole real hne. The periodic density profiles discussed in this paper 
are certainly not of this type. Nevertheless, they arise correctly as solutions of the 
variational equation associated with the alternative form of the collective potential 
containing the p^{x) term, given in Eq.(B.2), as they do, for example, in the pioneer- 
ing work [15], where these periodic density waves were discovered. The resolution 
of this paradox lies in proper readjustment of the chemical potential enforcing the 
constraint (7). 

Finally, in order for this paper to be self-contained, we briefly summarize in 
Appendix C the derivation of the collective- field Hamiltonian (2) from (1). 

2 Periodic BPS Density Waves: Soliton Crystals 

The Hamiltonian (2) is essentially the sum of two positive terms^. Its zero-energy 

classical solutions are zoro-momentum, and therefore time independent configura- 

^Recall the constraint (7) and our comment concerning H^ing following (7). In addition, as was 
mentioned above, the external confining potential was set to zero. Thus, the first two terms in (2) 
comprise the BPS limit of the model. 
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tions of the collective field (4), which are also solutions of the BPS equation 

S[p] = ^^ + A-/^ = 0. (16) 
2 p J X — y 

It is easy to check that the duality transformation (8) maps a solution p(x) of 
(16) with coupling A onto another solution p{x) — —\p{x) of that equation with 
coupling A = ^. As we shall see below in Eq. (33), all solutions of (16) are of definite 
sign, and never vanish along the real axis. Thus, such a positive solution of (16) is 
mapped by (8) onto a negative solution, and vice-versa. 

The BPS equation (16) may be written alternatively as 

(A-l)a,p = 27rApp^, (17) 

where p^ is the Hilbert-transform (A.l) of p. Note that for A = 1, where the CM 
describes non-interacting fermions^ the only solution of (17) is p = po = const. 
Henceforth, we shall assume A 7^ 1. 

The proper way to solve this nonlinear integro-differential equation is to consider 
it together with its Hilbert-transform[18, 28] 

(A - 1) d^p" = 7rA((p^)^ - + pI) , (18) 

where on the RHS we used the identity (A. 14) (and the fact that d^p^ = {dxp)^ on 

the LHS) . Here po is a real parameter such that 

00 

I dx{p{x)-po)^0. (19) 

—00 

It arises from the fact that we seek a solution of p{x) which need not necessarily 
decay at spatial infinity. (See (A. 2 ).) Note that (18) is even in po. By definition. 



^The constant solution is also the sole solution of (17) when A = 0, corresponding to non- 
interacting bosons. 
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the sign of po coincides with that of p{x), the solution of (18). A positive solution 
p{x) > corresponds to a BPS configuration of particles, and a negative one - to a 
configuration of antiparticles, as was mentioned following (8). 

To arrive at the Riccati equation mentioned in the introduction, we proceed as 
follows. Given the density p{x), consider the resolvent 

— oo 

associated with it, in which z is a complex variable. It is easy to see that it is related 
to the resolvent G{z) of the subtracted density p{x) = p{x) — po, defined in (A. 3), by 

^z) = - - G{z) + ipo sign i^z) . (21) 

The resolvent ^{z) is evidently analytic in the complex plane, save for a cut 
along the support of p{x) on the real axis. From the identity (A. 4) we obtain 

^±{x) = ^x ± iO) = p" {x) ± ip{x) , (22) 

consistent with (21) and (A. 5). Thus, if $(2;) is known, p{x) can be determined from 
the discontinuity of ^{z) across the real axis. 

An important property of $(2;), which follows directly from the definition (20), 

is 

^>$(.) = — /f^. (23) 
TT J \z — yr 
-00 ' 

Thus, if p{x) does not flip its sign throughout its support, we have 

sign (Q^<l>(z)) = sign (S>2;)sign (p(x)) . (24) 

We shall use this property to impose certain further conditions on the solution of 
(26) below. 
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It follows from (22) that (17) and (18) are, respectively, the imaginary and real 
parts of the Riccati equation 

(A - l)d,^±{x) = 7rA($|(x) + pI) (25) 

obeyed by both complex functions ^±{x) . Let ^±{z) be the analytic continuations 
of ^±{x) into the z— upper and lower half planes, respectively. These functions are 
evidently the two solutions of 

{x-i)dM^) = ^Km' + pl), (26) 

subjected to the boundary conditions ^^{x+iO) = $_(;r— iO) and sign ($5$+(.t + iO)) 
sign {p{x)) — signpo, from (22). The resolvent (20) is then obtained by patching 
together ^+{z) in the upper half-pane and ^-{z) in the lower half-plane. 

The standard way to solve (26) is to write it as 
f I 1 



dMz)-^k, (27) 



where 

h = 

A-1 ' 

is a real parameter. 



k=P^. (28) 



Straightforward integration of (27) then yields the solutions 

1 _)_ pikz-u± 

^^{z) = zpo ^ _ , (29) 

where u± are integration constants. The boundary condition $^(a; -|- iO) = $_(x — 
iO) then tells us that U- — —u*^. Clearly, Qu+ can be absorbed by a shift in 
x. Therefore, with no loss of generality we set — 0. The second boundary 
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condition sign (S5<l>+(x + iO)) = sign po then tells us that u = > 0. Thus, 

$±(z) are completely determined and we obtain (20) as 

\ _j_ gifez— usign (Qfz) 
*(^) = ^_g,fe,-„signO.) ■ (30) 

As can be seen in (33) below, the density p{x) associated with (30) is indeed of 
definite sign, namely, signpo. 

The asymptotic behavior of (30) is such that 

$(±ioo) = iiposign/c . (31) 

This must be consistent with (24), which implies (together with the fact that 
sign {p{x)) = signpo) that k must be positive. In other words, as can be seen 
from (28), positive (space-dependent) BPS density configurations (po > 0) exist 
only for A > 1, and negative (space-dependent) BPS densities (po < 0) arise only for 
< A < 1"^ . The duality symmetry (8), which interchanges the domains < A < 1 
and A > 1, maps these two types of BPS configurations onto each other. 

Now that we have determined ^{z), let us extract from it the BPS density p{x) 

and its Hilbert transform p^{x). Prom (30) we find that 

^ , , ^, „, — sin kx + isinhu 

X = $ X + iO) = Po r— , 32 

cosh u — cos kx 

from which we immediately read-off the solution of the BPS-equation (16) as 

, , sinhu 
P{x) = Po 



cosh u — cos kx 



H/ \ sin kx /oo\ 

P ^ P^ cosh u — cos kx ' 



^Constant solutions p = po of (16), are of course not subjected to this correlation between 
sign Po and the range of A. 
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where both A; > and u > 0, and the sign of p{x) coincides with that of po- That 
in (33) is indeed the Hilbert-transform of p can be verified by exphcit calculation. 

The static BPS density- wave, given by p{x) in (33), is nothing but the finite- 
amplitude solution of [15]. It comprises a two-parameter family of spatially periodic 
solutions, all of which have zero energy density, by construction. The period is 

^ 27r A-1 

T= — = -— . 34 

k Xpo 

It can be checked by explicit calculation^ that 

^ J p{x) dx = Po, (35) 

period 

oo 

and therefore that / {p{x) — po)dx — , as required by definition of po. Thus, 

—oo 

the parameter po determines both the period of the solution p(x), as well as its 
period-average, and the other (positive) parameter u determines the amplitude of 
oscillations about the average value. Note also from (35), that the number of parti- 
cles per period is 

Tpo = ^ . (36) 

A couple of limiting cases of (33) are worth mentioning. Thus, if we let m ^ , 
we obtain a comb of Dirac 5— functions 

p{x) = ^ E -^(^ - nT) . (37) 

If, in addition to m ^ 0, we also let k tend to zero (or equivalently, let the period T 
diverge), such that & = f remains finite, we obtain the BPS soliton solution [15, 16] 

/ N A- 1 1 b 



^The best way to do this computation is to change variables to t = e'*^^ and transform the 
integral into a contour integral around the unit circle. 
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In fact, the original construction of the periodic sohton (33) in [15] was done by 
juxtaposing infinite sohtons hke (38) in a periodic array. For this reason we may 
refer to the finite amphtude BPS density wave in (33) also as the soliton crystal. 

Note that the relation (36) is preserved in both limiting cases, since the RHS of 
(36) depends neither on u nor on k. 



2.1 BPS Solutions of the Two-Family Model at the Special 
Point (14) 



The BPS-equations of the two- family collective field Hamiltonian (12) are 



Solutions of these coupled equations yield the time-independent zero-energy and 
zero- momentum configurations of the collective fields pi and p2- 

Finding the general solution of these coupled equations for arbitrary couplings 
and masses (subjected to (10)) is still an open problem, which we do not address in 
the present paper. However, at the special point (14), where the two-family model 
becomes similar to a single-family model, the two equations (39) simplify drastically, 
becoming linearly dependent. For example, for 





(39) 



A — Ai — -— , Ai2 



1, 



and 



m — mi — — Am2 



(40) 



it is easy to see that 




(41) 
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Since at the same time, from (39), Bi = B2 = 0, (41) implies that the two densities 
must be proportional 

P2{x) = -Kpi(x) . (42) 

(From the discussion in [4] we know that the constant k > 0, and the negative 
density is interpreted as density of holes, as was mentioned in the Introduction.) 
Upon substitution of (42) back in (39) we see that 

Bi = — Xirp.ff , (43) 

^ Peff 

coincides with the corresponding one-family expression B in (16) with an effective 
density pe// given by (15). Thus, at this special point, peff{x) is given by (33), from 
which pi and p2, being proportional to Peff, can be deduced as well. An analogous 
solution of (39) exists for the case in which the roles of the two families in (40) are 
interchanged. 

3 Non-BPS Solutions of the Equation of Motion 

The uniform-density ground state, as well as the periodic space-dependent BPS- 
configurations discussed in the previous section, all correspond to zero-energy and 
zero- momentum configurations of the collective field Hamiltonian (2) . Static density 
configurations with positive energy density are found by extremizing the collective 
potential 

- ^jdx p{x)B[pf + ,,{n- j dx p{x)) (44) 
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part of (2). 

Computation of the variation of (44) with respect to p is most easily carried with 
the help of (45) just below. Thus, using the elementary relation / dx F{x) G^{x) = 
— J dx F^{x) G{x) it is easy to obtain the variational identity 

A- 1 



j dx p{x) F{x) 5B[p] = j dx 
where the infinitesimal variation of B[p\ 



2p 



d,{pF) + 7iX{pF) 



H 



6p{x) , (45) 



5B[p]^^d,(^^^ -nXSp"" (46) 

was computed from (16). 

Using (45), it is straightforward to obtain the desired variational equation as 

2m B[p]' - ^ d,{pB[p]) + 2n\ {pB[p]f - 2m// = . (47) 

opyx) p 

The collective potential (44) is invariant under the duality transformation (8). 
Thus, (47) must transform covariantly under (8). Indeed, it is straightforward to see 
that under (8), the variational equation (47) transforms into times itself. In this 
way, a solution p{x) of (47) with parameters A, m, ji is transformed into a solution 
p{x) of (47) with parameters A, m, jl. We shall make use of this fact later-on. 

Evidently, any solution of the BPS equation B[p] — (Eq.(16)) is also a solution 

of the variational equation (47) (with = 0), reflecting the fact that (44) is quadratic 
and homogeneous in B[p\. 

Unfortunately, we do not know how to find the most general solution of this 
equation. Therefore, we shall content ourselves with finding a particular family of 
solutions to (47) of the simple shifted form 

p{x)^ Ps{x)^c, (48) 
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where 



, , sinhtt 
Ps\X) cosh, u — cos kx 



is the BPS profile in (33) and c is an unknown constant, to be determined from (47). 
(Clearly, (48) with c = must be a solution of (47).) 

Let us proceed in a few steps. First, note that 

from (33). Then, compute 

\-ld^p H Xncposinkx 

B[p\ = — Xnp — —. ; r — — , (51) 

2 p c(cosh u — cos kx) + po smh u 

from which we obtain the remarkably simple relation 

pB[p\ = -Xttcp^ . (52) 

Therefore, 

{pB[p]f ^Xncip.-po), (53) 

where we used the identity (A. 8). (Note that (53) is consistent with the identity 
F^{x) dx — 0.) Substituting the ansatz (48) and the auxiliary results (51)- 
(53) in (47) we obtain the LHS of the latter as a rational function of polynomials in 
cos kx. The numerator of that function is a cubic polynomial, which we then expand 
into a finite cosine Fourier series, all coefficients of which must vanish. Thus, the 
coefficient of cos ?>kx determines the chemical potential in terms of the remaining 
parameters as 

/^ = -^^Po(po + 2c), (54) 
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which we then feed into the coefficients of the remaining three terms. The coefficient 
of cos2A;a; is then found as the cubic 

— (A7r)^Po c(c^ sinh u + 2cpo cosh u + pg sinh u) , (55) 

where we used (28) on the way. This coefficient must vanish, yielding a cubic equa- 
tion for c. The remaining Fourier coefficients vanish identically upon substitution 
of the roots of this cubic equation for c. 

As we have anticipated following (48) , one root of this cubic equation is obviously 
Co = 0, which corresponds to p = p^. The other two roots are 

Ci = — po tanh — and C2 = — po coth — . (56) 

Note that neither of these roots, and therefore neither of the shifted solutions (48), 
depend on m or on /i. Once the parameters m and fj, are related according to (54), 
they drop out of any further consideration. 

3.1 Large Amplitude Density Waves: Vortex Crystal Solu- 
tions 

From this point onward we shall discuss the cases A > 1 and < A < 1 separately. 
3.1.1 The case A > 1 

In this case po > 0, as we saw following (31). For positive po, the ffist root Ci in (56) 
amounts in (48) to shifting the BPS solution Ps(x) by its minimum. The resulting 
solution 

/ \ f sinhw u\ . 

P^^'^^P' Uh^-cosA:. -^""^ 2 J 

20 



is a positive function which vanishes periodically. We shall refer to it as the vortex 
crystal solution, as it is a periodic generalization of the single vortex solution of 
[16]. Since Pp{x) > , it is a density of particles (rather than holes). Therefore it 
corresponds to having a positive mass parameter m > in (44). The vortex crystal 
(57) coincides with the so-called large amplitude wave solution of [17] for the case 
A > 1 and zero velocity. 

The second root C2 in (56) amounts to shifting the BPS solution Ps{x) by its 
maximum. The resulting solution 



is thus a negative function which vanishes periodically - an anti- vortex crystal. Being 
a negative solution of (47), Eq.(58) should be interpreted as the density of holes 
rather than particles. Therefore it corresponds to having a negative mass m < in 



In this case po < 0, as we saw following (31). Therefore Ci and C2 in (56) switch 
roles: For negative po, ci amounts to shifting in (48) the BPS solution Ps{x) by its 
maximum. The resulting solution 



is a negative function which vanishes periodically - an anti-vortex crystal. It is 
therefore a density of holes corresponding to having a negative mass m < in (44). 

The second root C2 amounts in this case to shifting the BPS solution Ps{x) by 




(58) 



(44). 



3.1.2 The case < A < 1 




(59) 
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its minimum. The resulting solution 
_ / X / sinhti 




— coth — = IpoI coth 



sinli u 



(60) 



cosh u — cos kx 



cosh u — cos kx 



is thus a positive function which vanishes periodically - a vortex crystal. It cor- 
responds to having m > in (44), in a similar manner to Pp{x) in (57). Pp{x) 
coincides with the large amplitude wave solution of [17] for the case < A < 1 and 
zero velocity. Note that Pp{x) has appeared also in [18]. 

The duality transformations (8) leave the wave-number k in (28) invariant. By 
definition, the positive parameter u, defined in (30), is invariant under (8) as well. 
Thus, evidently, the duality transformations (8) map Pp{x) in (57) and Pn{x) in (59) 
onto each other. (Of course, the po parameters appearing in the latter two equations 
are different from each other, and related by the fourth relation in (8).) Similarly, 
the duality transformations (8) map Pn{x) in (58) and Pp{x) in (60) onto each other. 

3.1.3 Average Energy Densities per Period 

Our new solutions (57) - (60) of the variational equation (47) are periodic functions, 
with the same period as that of the BPS solution pg. Since these are non-BPS 
configuration, they must carry positive energy density®. We shall now proceed to 
calculate the mean energy densities per period of these configurations, to which end 
we must determine the combination pB[p\^ appearing in (44). Prom the general 
expressions (51) and (52) we obtain 



^Since negative densities eorrespond to holes, whose mass should be taken negative, we have 
> in these cases as well. This renders VcoU in (44) positive for such densities. Thus, the 
negative solutions p„ and p„ carry positive energy density, as their positive counterparts obviously 
do. 



pBip? 



Xncpo sin kx Avrcpo sin kx 



(61) 



cosh u — cos kx c(cosh u — cos kx) + po sinh u 
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The desired period-averaged energy density is then given by 



£ 



dx 



period 



2m 



(62) 



where T- f (Eq. (34)). 



We shall content ourselves with computing the energy density only of the positive 
densities Pp and Pp. In order to compute S of (57), corresponding to A > 1 and 
Po > 0, we substitute c = Ci in (61). After some algebra, we find that in this case 



u 



PpB[pp\ = (Attpo) tanh - 



u zt 
Po(l - tanh -) - {ps{x) - po) tanh - 



(63) 



In view of (35), and by definition of po, the period-average of Ps{x) — po vanishes. 
Thus, from (44), we obtain the period-average energy density of (57) as 



(Attpo) 



u 



Po tanh —(1 — tanh — ) 
2m ^ 2^ 2' 



(64) 



which is a manifestly positive quantity. It depends continuously on the two param- 
eters Po a-nd u, comprising an unbounded continuum of positive energies, which is 
not gapped from the zero energy density of the BPS sohtons. 

Similarly, in order to compute £ of (60), corresponding to < A < 1 and po < 0, 
we substitute c = C2 in (61). After some algebra, we find that in this case 



PpB[pp] = -(Attpo) coth 



u 



Po(coth I - 1) + {ps{x) - Po) coth I 



(65) 



which leads to the positive period-average energy density of (60) given by 



^[Pp] = - 



(Avrpo) 



u 



u 



2m PO coth - (coth --1) 



(66) 
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3.1.4 Energy Densities at Fixed Average Particle Density 



It is particularly useful to consider the energy densities (64) and (66) at a fixed 
average particle density per period. The latter is, of course, the subtraction constant 
as defined in (19), which is given by 

po = Po + c (67) 

for the shifted solutions (48). Both (64) and (66) depend on the two parameters 
Po and u. Holding po fixed can thus be used to eliminate one of these parameters, 
which we shall take to be u. 

Let us concentrate first on pp in (57), for which c — ci — —po tanh | (and of 
course, A > 1). Thus, 

u 

Po = Po (1 - tanh -) , (68) 

which is positive, since po > in (57). Moreover, po > po in this case, since u > 0. 
In terms of this fixed po, we obtain S[pp] in (64) as 

^[Pp] = PoPo (po - Po) , Po > Po = fixed . (69) 

2m 

This energy density vanishes at the minimal possible value of po = po, corresponding 
to « = 0, and therefore to the BPS density configuration (37). As po increases from 
its minimal value, the period-average energy density £[pp] increases monotonically 
from zero to infinity. Increasing po really means increasing the wave number k — 
^^3^, i.e., making the density modulation wave-length shorter. 

It is interesting to note that in terms of k and po we can write 

1 (l-A)po fl-X 2 , Avrpo \ 27rA 27rA . 

£[pp] = k + -^kj , k = ^po > j—^po = fixed , (70) 
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where the expression within the brackets is nothing but the dispersion relation for 
fluctuations around the constant background po [12]. 

We can analyze the periodic vortices pp in (60) in a similar manner. For these 
solutions C — C2 — —P0 coth | (and of course < A < 1). Thus, 

Po = Po (1 - coth -) = IpoI (coth 2 ~ 1) ' C^^) 

which is again positive, since the allowed range of po in (60) is po < 0. For a given 
value of Po, po = ~^{^^ ~ l)Po ranges throughout the negative real axis as u ranges 
throughout the positive one. In terms of this fixed po, we obtain an expression 
for S[pp] in (66) which coincides with the RHS of (69), but where now po < 0, of 
course: £[pp\ — Pol Pol (I Pol + Po) ■ This energy density vanishes at the maximal 
possible value of po = 0, corresponding to u — 0, and therefore to the BPS density 
configuration (38). As po becomes increasingly negative the period-average energy 
density £[pp\ increases monotonically from zero to infinity. In terms of the wave 
number k and po, we obtain that S[pp\ = k'^ + ^^A;) , which coincides 

with (70), but where now A; > for any value of po > 0. 

3.2 The Two-Family model at the Special Point (14) 

The variational equations associated with the two-family collective potential part of 
the two-family collective Hamiltonian (12) are 

Bl -h^d^^p^B^)^ 2X^7: {p^B^f + 2^X^27^ {p^B^f -2m^p^ ^ 
Pi ^2 

2_A2^^^ 2A27r(p2S2)^ + 2— Ai27r(pi5i)^-2m2P2 = 0, 

P2 rni 

(72) 
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in straightforward analogy with (47), where the BPS combinations Bi and B2 were 
defined in (39). 

As in the case of the BPS equations (39), the general solution of these coupled 
equations for arbitrary couplings and masses (subjected to (10)) is still an open 
problem, which we do not address in the present paper. However, at the special 
point (14), where the two-family model becomes similar to a single-family model, 
the two equations (72) simplify drastically, becoming linearly dependent, in much 
the same way that the BPS equations (39) got simplified. 

Consider, for example, the case (40). In this case (41) still holds, of course, 
but now neither Bi nor B2 need vanish. Thus, we cannot conclude that pi and p2 
must be proportional. Instead, we shall now show that under the condition (40), 
there is a non-BPS solution of the coupled equations (72) in which the two densities 
are proportional to each other, as in (42). In this case it follows from (41) that 
Bi + XB2 = . Substituting this relation as well as (40) in (72) , we see that the two 
equations coincide, provided pi + XiJ,2 — 0, and that their common form is nothing 
but the variational equation (47) of the single-family, for an effective density (15). 
Thus, at this special point, Peff{x) is given by (49), (57) or (58), from which pi and 
P2, being proportional to Pe// can be deduced as well. An analogous solution of (72) 
exists for the case in which the roles of the two families in (40) are interchanged. 
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Appendix A: A Compendium of Useful Hilbert- 

Transform Identities 



For the sake of completeness, and also for future reference, in this Appendix we list 
and prove some well-known and useful identities involving Hilbert transforms. 

Consider the class of (possibly complex) functions p{x) on the whole real line 
— oo < X < oo, whose Hilbert transforms 

= If iy£^ (A.1) 

exist, and which can be made integrable by subtracting a constant po- Let us denote 

p{x) = p{x) - po . (A.2) 

(If p{x) is already integrable, then po = 0, of course.) Thus, for example, if 

oo 

p{x) is periodic with period T, with a Fourier zero-mode po, then / dxp{x) — 

— oo 

oo 

/ dx {p{x) - Po) = . 

— oo 

Given p(x), consider the resolvent 

G{z)= J (A.3) 



-oo 



associated with it, in which 2; is a complex variable. The resolvent G{z) is evidently 
analytic in the complex plane, save for a cut along the support of p{x) on the real 
axis. From the identity 

^ ^ ± i7rS{x) , (A.4) 



x^iO X 

where P denotes the Cauchy principal value, we then obtain the well-known formula 

G{x T iO) = -Trp"{x) ± iTTpix) , (A.5) 
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where in the term one before last we used the fact that p^{x) = p^{x). Thus, if 
G{z) is known, p{x) can be determined from the discontinuity of G{z) across the 
real axis. 

As a nontrivial example consider p{x) — p{x) — e*^. For this function p^{x) — 
ie^^ = ip{x) and G{z) = -2m e{Qz) Consequently G{x-iO) ^ and G{x+iO) = 
— 27rie*^, in accordance with (A. 5). As yet another example consider the Cauchy 
probability distribution p{x) — p{x) — ^p^^. For this function p^{x) — —^ ^i^^a 
and G{z) = ^_|_;^sign (Qz) • Consequently G(x^iO) = in accordance with (A. 5). 

For all functions in this class, as 2; — >■ oo, G{z) tends asymptotically to zero not 
slower than -, that is 

z ' 

G(z) ^ o(-) . (A.6) 

00 

If, in addition, all moments M„ = / dxx" p{x) , (n > 0) of p{x) exist, then G{z) 

—00 

is the moment generating function of p{x), namely, it has the large- ;s expansion 




The analyticity properties of G{z) and the bounds on its asymptotic behavior 
at infinity are at the heart of our derivation of the Hilbert-transform identities to 
follow. 

From this point on, we shall take all functions p{x) to he real. For real p{x), we 
deduce from (A. 5) that 

pH(^x) = --UG{x T ^0) and p{x) = ±-^G{x ^ «0) . (A.7) 

TT TT 

As a warm-up exercise, let us prove the well-known fact that 

{p''{x)f ^--p{x)^p^-p{x). (A.8) 
28 



Thus, consider 




'^{y - ) {p^{y) + ip{y)) dy - p{x) 



where in the last step we used (A. 4) and (A. 7). Let us now prove that the last 
integral vanishes, from which (A. 8) follows. To this end, complete the contour of 
integration in the last integral (namely, the line running parallel to the real axis 
just above it) by the infinite semi-circle in the upper half-plane > 0, traversed 
in the positive sense. Let us denote the closed contour thus formed by 7. Due to 
the asymptotic behavior (A. 6) of G{z) we can establish the first equality in 



whereas the second equality follows since the contour 7 encompasses no singularity. 
We shall now prove the important identity 



obeyed by any two functions pi{x) and P2{x) in the class of functions considered. 
Our first step in proving (A. 9) is to observe that it may be written equivalently as 




—00 



{P1P2 + Pi P2) = Pi P2 - P1P2 + P10P20 



(A.9) 



(pipf + PiP2r = pf pf - P1P2 . 



(A.IO) 



Consider now the contour integral 




(A.ll) 
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where Gk{z) is the resolvent corresponding to pk{x) {k = 1,2) , x E M, and where 
Coo is the circle of infinite radius, centered at the origin. Due to the asymptotic 
behavior (A. 6) of the two resolvents, evidently 

7 = 0. (A.12) 

Since Gi^2{z) are analytic off the real axis, we can deform Coo into the positively 
oriented boundary F of an infinitesimal strip around the real axis (namely, the union 
of a line parallel to the real axis just below it, traversed from — oo to oo, with a line 
parallel to the real axis just above it and traversed in the opposite direction). The 
contour integral around F essentially picks up the imaginary part of the integrand 
evaluated just above the real axis. Thus, we have 



/r z — X 2711 TT J y — X + iO 

— oo 

The last integrand may be written as 

7T^^^-iiT5{y-x)] n (Pkiy) + miy)), 

by virtue of (A. 4) and (A. 7). Upon substituting the last expression in (A.13) and 
taking the imaginary part, we obtain the desired result (A. 10). Note that for pi = 
P2 — p, (A. 9) simplifies into 

2ipp^f = {p^r-p' + pl (A.14) 

Finally, we shall prove an identity involving three functions Pk{x) {k — 1,2,3) 
and their Hilbert transforms. Our proof follows essentially the one given in [29, 30] 
for the case pi = p2 = P3, which is reproduced also in the text-book [8]. Let Gk{z) 
be the resolvent corresponding to Pk{x). Consider now the contour integral 

J=f ^Gi{z)G2{z)Gs{z). (A.15) 
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As in the previous proof, due to the asymptotic behavior (A. 6) of the resolvents, 
evidently 

J = 0. (A.16) 

Since the resolvents are analytic off the real axis, we can deform Coo into the contour 
r, as in the previous proof, which picks up the imaginary part of the integrand 
evaluated just above the real axis. Thus, we have 

1 °° 

= J= — » [ Gi{y + iO)G2{y + iO)G3{y + tO) dy . (A.17) 

— oo 

The last integrand may be written as 

II {pk iy) + m{y)) , 

k=l 

by virtue of (A. 7). Upon substituting the last expression in (A.17) and taking the 
imaginary part, we obtain the desired result 

oo oo 

/ {p?P^P3 + P?P2P3 +P1P2P3) J PiP2P3dx. (A.18) 

—00 —00 

Note that for pi — p2 — p^ — p, (A.18) simphfies into 

00 00 
2, j pip'^fdx^ j {pfdx, (A.19) 
—00 —00 

which is the identity proved in [8, 29, 30] . 

Since (A.18) holds for any triplet of functions pk in the class of functions thus 
considered, we can write it formally as an identity among distributions acting upon 
these test functions, namely, the well-known [31] identity 

P P P P P P 

— — + — — + — — = TT^dix - y)d(x - z) . (A.20) 

x—yx—z y—xy—z z—xz—y 
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In [31], the identity (A. 20) was proved using Fourier transforms. For alternative 
proofs of the identities discussed in this Appendix, and for more information about 
Hilbert-transform techniques, see Appendix A of [27]. 

As should be clear from our proof, (A. 20) holds only when the distributions on 
both its sides act upon functions which are integrable on the whole real line. How- 
ever, this identity is frequently used in the literature on collective field theory beyond 
its formal domain of validity. For further discussion of this problem see Appendix 
B, where we show that this transgression is benign, and can be compensated for by 
readjusting the chemical potential which governs the normalization condition (7). 
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Appendix B: A Paradox and its Resolution 



The expression for the collective potential in (44) contains bilocal as well as trilocal 
terms in the density. It is customary in the literature to avoid the trilocal terms by 
applying a standard procedure as follows: The principal value distribution, acting 
on functions integrable along the whole real line, satisfies the identity (A. 20), which 
we rewrite here for convenience 

P P P P P P 

+ — — + — — = 7r^S{x - y)5{x - z) . (B.l) 



x—yx—z y—zy—x z—xz—y 
Making use of (B.l) in (44), we obtain 

i^-r f...s , (A-i)' A.JM' , MA-i) f,p{y) 



Vcoll — 



bm J am J p 2m J J x — y 



+ p{n- j dxp(x)j (B.2) 

This expression for Vcou is evidently devoid of any trilocal terms. (Note that the 
chemical potential fl in (B.2) need not coincide with the one in (44), as our notations 
imply) 

The classical equation of motion which results from varying (B.2) is 



2m 8m p Am p m J x — y 

It was this form of the equation of motion (rather than (47)), from which the solitons 
and density waves were derived in the pioneering work [15]. 

It can be checked that ps in (49), pp in (57) and p„ in (58), the solutions of the 
variational equation (47) of the first form (44) of the collective potential, are also 
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solutions of (B.3) (albeit, with values of jl different from those of (54)). That this 
is true may look surprising, and even paradoxical to some readers, since neither of 
these solutions is integrable along the whole real line, which is a necessary condition 
for (B.l) to hold. This should be clear from the proof of (A. 20) in Appendix A, but 
it can also be demonstrated by a simple counter example - just apply both sides 
of (B.l) on three constant functions and integrate over all coordinates. The LHS 
would be null, while the RHS would diverge. 

In fact, the latter counter example is precisely relevant to determining the ground 
state of the collective Hamiltonian (2). The uniform ground state density p = Qq is 
a solution of the BPS equation (16), and of course, also of the variational equation 
(47) with fjL — O. The energy density tied in it is of course null. It is also a solution 
of the alternative variational equation (B.3) with jl = '■'^^'^"^ and energy density 



Thus, it seems that using (B.l) beyond its formal domain of validity is a mild 
transgression, which is compensated for by appropriately readjusting the chemical 
potential. This is indeed true, as we shall now prove, thus resolving the paradox why 
(47) and (B.3) always lead to the same solutions. To this end we shall consider all 
p configurations which are simultaneous solutions of (47) and (B.3). Such functions 
are evidently extrema of AV = Vcoii — Vcoii- From (44) and (B.2) we obtain 



(Note that we have omitted from this expression the constant term {p — p)N .) Due 
to the analytic structure of $(-2), and as explained in Appendix A, the latter integral 



(with respect to (B.2)) 



6m 





(B.4) 
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can be written as the contour integral 

where is a circle of infinite radius, centered at the origin. Note that ^{z) need 
not decay as 2; ^ 00, since J dx p may diverge. Thus, in general AV ^ 0. 



We shall now determine solutions of 

SAV 



Sp{x) 

To this end, let us compute 

6^{z) 15 J p{u) du 



= 0. (B.6) 



5p{x) TT Sp(x) J z — u n z — X 
—00 

Prom this we infer that 

6AV (A7r)2 r dz ^^{z) 



(B.7) 



op[x) 2m Jcoo 27^^ z — x 
The contour Coo in the last integral can be deformed to the countour F, defined in 
Appendix A, which essentially picks up the imaginary part of the integrand evaluated 
just above the real axis. Thus, in a manner similar to the discussion in Appendix 
A, from (A.9) to (A. 14), we obtain 

But from the identity (A. 14) we see that the latter equation boils down to 

5AV _ {Xnpof 

In other words, the condition (B.6) simply relates the two chemical potentials 

(B.n) 
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setting no further conditions on p{x), where po is the subtraction constant associated 
with the p in question, and should not be confused with the one appearing in (49). 

To summarize - any solution of (47) with chemical potential /j, is simultaneously 
a solutions of (B.3) with chemical potential jl given by (B.ll). 
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Appendix C: A Brief Summary of the Collective 

Field Formulation of the Calogero 
Model 



In order for this paper to be self-contained, we briefly summarize in this appendix 
the derivation of the collective- field Hamiltonian (2) from (1) . 

The singularities of the Calogero-model Hamiltonian (1), namely, 

at points where particles coincide, implies that the many-body eigenfunctions con- 
tain a Jastrow-type prefactor 

TV 

n = n (^^ - . (C.2) 

This Jastrow factor vanishes (for positive A) at particle coincidence points, and 
multiplies that part of the wave-function which is totally symmetric under any per- 
mutation of particles^. It is precisely these symmetric wave- functions on which the 
collective field operators act, as explained below. 

Let us recall at this point some of the basic ideas of the collective-field method 
[8, 9, 10], adapted specifically to the Calogero model[ll, 12]: Instead of solving 
the Schrodinger equation associated with (C.l) for the many-body eigenfunctions, 
subjected to the appropriate particle statistics (Bosonic, Fermionic of fractional), we 
restrict ourselves to functions which are totally symmetric under any permutation 
of identical particles. This we achieve by stripping off the Jastrow factor (C.2) from 



''Note, in particular, that for A = and A = 1, the model describes interacting bosons and 
fermions, respectively. 
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the eigenfunctions, which means performing on (C.l) the similarity transformation 

H ^ H = Il-^HIl, (C.3) 

where the Hamiltonian 

1 ^ r)2 \ N 1 o 

2m ^ dxj^ m ^ - Xj dxi 

Note that if does not contain the singular two-body interactions. By construction, 
this Hamiltonian is hermitian with respect to the measure 

dix{xi) = d^x , 

(as opposed to the original Hamiltonian H in (C.l), which is hermitian with respect 
to the flat Cartesian measure) . 

We can think of the symmetric many-body wave-functions acted upon by if as 
functions depending on all possible symmetric combinations of particle coordinates. 
These combinations form an overcomplete set of variables. However, as explained 
below, in the continuum limit, redundancy of these symmetric variables has a neg- 
hgible effect. The set of these symmetric variables can be generated, for example, 
by pro dues of moments of the collective - or density - field 

N 

p{x) = '^5{x - Xi) . (C.5) 

i=l 

The collective-field theory for the Calogero model is obtained by changing variables 
from the particle coordinates Xi to the density field p{x) . This transformation 
replaces the finitely many variables Xi by a continuous field, which is just another 
manifestation of overcompleteness of the collective variables. Clearly, description 
of the particle systems in terms of continuous fields becomes an effectively good 
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description in the high density hmit. Of course, the large density hmit means that 
we have taken the large- N limit. 



Changing variables from particle coordinates Xi to the collective fields (C.5) 
implies that we should express all partial derivatives in the Hamiltonian H in (C.4) 
as 

dxi I dxi 5p{x) ' ^ ^ 

where we applied the differentiation chain rule. 

In the large — A'^ limit, the Hamiltonian H can be expressed entirely in terms of 
the collective field p(x) and its canonical conjugate momentum 

.M = -,^. (C.7) 

as we show below. It follows from (C.6) and (C.7) that the particle momentum 
operators (acting on symmetric wave-functions) may be expressed in terms of the 
collective-field momenta at particular points on the line as 

Vi = (C.8) 

(where 7r'(x) = dxT^{x)). Finally, note from (C.5) that the collective field obeys the 
normalization condition 

j dxp{x) = N . (C.9) 

The density field p and its conjugate momentum tt satisfy the equal-time canonical 
commutation relations® 

[p{x),TT{y)]^i5{x-y), (C.IO) 



^According to (C.9), the zero-momentum modes of the density fields are constrained, i.e., non- 
dynamical. This affects the commutation relation (C.IO), whose precise form is [p{x),'K{y)] = 
i{5{x — y) — (1/0)) where I is the size of the large one-dimensional box in which the system is 
quantized, which is much larger than the macroscopic size L of the particle condensate in the 
system. In what follows, we can safely ignore this l/l correction in the commutation relations. 
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(and of course [p{x),p{y)] = [7r{x),7r{y)] = 0) . By substituting (C.5)-(C.8) in (C.4), 
we obtain the continuum-limit expression for H as 

(C.ll) 

where f- denotes Cauchy's principal value. 

It can be shown [8] that (C.ll) is hermitian with respect to the functional 

measure^ 

Vpi[p]^J[p]l[dp{x), (C.12) 

X 

where J[p] is the Jacobian of the transformation from the {xi} to the collective 
field {p{x)} . In the large - N limit it is given by [32] 

In J = (1 — A) J dxp{x) In p{x) ~ ^ J dxdyp{x) In \x — y\p{y) (C.13) 

It is more convenient to work with a Hamiltonian, which unlike (C.ll), is hermitian 
with respect to the fiat functional Cartesian measure dp{x) . This we achieve by 
means of the similarity transformation -0 J^ij) ,H Hcou = J^HJ~^ , where 
the continuum collective Hamiltonian is 

^ 2^ / ^'^^^ ^^^^ ^'^^^+2^ / ^^^^^^ + Vt^) 

(C.14) 

namely, the Hamiltonian given by (2) and (3). The collective-field Hamiltonian (12) 
of the two-family Calogero model can be derived from (9) in a similar manner. 
^By definition (recall (C.5)), this measure is defined only over positive values of p . 
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